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Abstract 

Diverse mass and mixing patterns between the quarks and leptons makes it challenging to con- 
struct a simple grand unified theory of flavor. We show that SO(10) SUSY GUTs with type II 
seesaw mechanism giving neutrino masses provide a natural framework for addressing this issue. 
A simple ansatz that the dominant Yukawa matrix (the 10-Higgs coupling to matter) has rank 
one, appears to simultaneously explain both the large lepton mixings as well as the observed quark 
flavor hierarchy in these models. A testable prediction of this ansatz is the neutrino mixing, U e ^, 
which should be observable in planned long baseline experiments. 
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I. INTRODUCTION 



Understanding the origin of the hierarchical pattern of quark masses and mixings has long 
been recognized as a challenge for physics beyond the standard model The discovery of 
neutrino masses and mixings with totally different flavor pattern than quarks (i.e. 8 l 23 ~ 45° 
and 8[ 2 — 35° as against 8 23 ~ 2.5° and 8\ 2 ~ 13°) has added more mystery to the flavor 
problem. In generic bottom-up pictures where quarks and leptons are treated as different 
species of particles with no particular relation between them, this problem is not so serious 
since one can simply focus on each sector separately, as is often done for neutrinos 
However, in grand unified theories where the quarks and leptons unify at a very high scale, 
one would naively expect that their masses and mixings would exhibit a similar pattern. The 
fact that they are so different may be hint of some really new exciting underlying physics. 
In this note we address this question in the context of supersymmetric SO(10) models with 
renormalizable Yukawa couplings being responsible for fermion masses. 

We show that in SO(10) models with 10, 126 plus possibly another 10 or 120 Higgs 
fields where fermion masses are g enerated by —lizable Yukawa coup.i„ g8 Q only and 
where type II seesaw is responsible for neutrino masses jffj], there is a natural way to have 
a unified understanding of both large lepton mixings and small quark ones. The basic idea 
is to require that one of the 10 Yukawa couplings is the dominant one contributing to up, 
down and charged lepton masses and has rank one with other smaller couplings providing 
neutrino masses as well as most of the quark lepton flavor hierarchy. Rank one plus small 
corrections as a way to unravel fermion flavor in D-brane models was discussed in jsj. We 
find that SO(10) models with type II seesaw are ideally suited for such an ansatz. 

A specific form of the rank one matrix can lead to tri-bimaximal mixing with corrections 
dictated by the quark flavor pattern. 

This paper is organized as follows: in sec. II, we review the mass formulae in SO(10) 
models with renormalizable couplings; in sec. Ill, we summarize our basic strategy for 
understanding the quark lepton flavor in a unified manner, discuss the rank one ansatz 
and apply it both the two generation case (IIIA) and three generation cases (IIIB). In sec. 
IIIC, we present realistic three generation models and outline their predictions. Sec. IV is 
devoted to some specific conjectures for the rank one matrix which could emerge from discrete 
symmetry models with specific discussion on the correction to tri-bimaximal mixings. Sec. 
V is devoted to a possible way to obtain the rank one ansatz in SO(10) models and in sec. 
VI, we present our conclusions. 
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II. OVERVIEW OF RENORMALIZABLE SUSY SO(10) MODELS FOR 
FERMION MASSES 



The basic idea in this class of models is to consider SUSY SO(10) theory with Higgs 



fields that give fermion masses to be in 10 (denoted by H) and 126 + 126 (denoted by 
A and A) plus either an extra 10 (H') or 120 (X). The GUT symmetry is broken by 
210+54+126+126 [8|. The Yukawa superpotential of this model is: 

W Y = h i/ji/jH + / ipipA + h! i/ji/j (E or H') , (1) 

where the symbol ip stands for the 16 dimensional representation of SO(10) that represents 
the matter fields. The coupling matrices h and / are symmetric, and h! is symmetric or 
anti-symmetric depending on whether we adopt H' or S. The representations H, H' and A 
have two standard model (SM) doublets in each of them whereas £ has four such doublets. 
The general way to understand so many SM doublets is that at the GUT scale Mjj, once the 
GUT and the B — L symmetry are broken, one linear combination of the up-type doublets 
and one of down-type ones remain almost massless whereas the remaining ones acquire GUT 
scale masses just like the color triplet and other non-MSSM multiplets. The electroweak 
symmetry is broken after the light MSSM doublets (to be called H u d ) acquire vacuum 
expectation values (vevs) and they then generate the fermion masses. The resulting mass 
formulae for different fermion masses are given by: 

Y u = h + r 2 f + r 3 ti, (2) 

Y d = n(h + f + ti), 

Y e = n (h-3f + c e ti), 
Y v d = h - 3r 2 f + c u ti, 

where Y a are mass matrices divided by the electro-weak vev v w k and and c e)V are the 
mixing parameters which relate the H u ^ to the doublets in the various GUT multiplets. 
More precisely, the matrices h, f and h! in Y a are multiplied by the Higgs mixings. The 
precise definitions of the couplings and the Higgs mixings are given in ref. 0]. When H' is 
adopted for the h! coupling, c e = 1 and c v = r 3 . In generic SO(10) models of this type, the 



neutrino mass formula has a type I [9| and a type II [10| contributions: 

M v = fv L -M D ^—M l Dl (3) 

J V R 

where vl is the vev of the B — L = 2 triplet in the 126 Higgs field and is given by vl ~ 
X ^r k ■ Note that in general, the two contributions to neutrino mass depend on two different 



parameters and it is easy to have symmetry breaking pattern in SO(10) ll| where the first 
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contribution (the type II term) dominates over the type I term. The neutrino mass formula 
then becomes 

M v = fv L . (4) 

Note that / is the same coupling matrix that appears in the charged fermion masses in 
Eq. (j2J), up to factors from the Higgs mixings and the Clebsch-Gordan coefficients. The 
equations (j2J) and (@J are the key equations in our unified approach to address the flavor 
problem. 

The main hypothesis of our approach is that the fermion mass formula of Eq. (j2J) are 
dominated by the matrix h with the contributions of / and h' being small perturbations. 
In the limit of /, h' — > 0, the quark and lepton mixings vanish as do the neutrino masses. 
We will show below that this simple hypothesis combined with Eq. fll]) can simultaneously 
explain large lepton mixings while keeping the quark mixings being proportional to 
and hence small. We will subsequently assume that the matrix h has rank one in which case 
the mass hierarchy can also be explained in a natural manner. 



III. EXPLAINING QUARK-LEPTON FLAVOR HIERARCHIES 

The quark and lepton mixing matrices are given by the product of diagonalizing uni- 
tary matrices for quark and lepton mass matrices as follows: denoting the diagonaliz- 
ing matrices of M u , Md by V u and Vd respectively (e.g., V u M u M]y^ = diag.(m^, ml, mf) 
and similarly for the down quark mass matrix ), the CKM (Cabibbo-Kobayashi-Maskawa) 
quark mixing matrix is given as V CKM = V U V\. The PMNS (Pontecorvo-Maki-Naka gawa- 
Sakata) lepton mixing matrix is given as £/pmns — (KKD* m the similar notation (e.g., 
V V M V V* = diag. (mi, m 2 ,m 3 )). 

In general, when two matrices with random 0(1) elements are considered, the mixing 
angles of the relative diagonalizing unitary matrices are all 0(1) in radian, while the eigen- 
values can have a hierarchy of 0(0. 1). In such an anarchical scenario, the neutrino masses 
and mixings can be explained (except for the CHOOZ bound of 13 neutrino mixing): the 
neutrino mixings are generically 0(1) and there is a little hierarchy for the neutrino mass 



squared difference ratio Am^/Am^ 12|. On the other hand, since the quark mixings are all 
smaller than 0(1) and the masses of quarks and charged leptons are very much hierarchical, 
anarchic mass matrices in general provide no explanation of these observations. Besides, it 
appears that the mass ratios and CKM mixings have several correlations among them. It is 
therefore to be expected that the quark and lepton matrices instead of being independent 
anarchic matrices must have some relations among them and an underlying theory leading 
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to these relations. In this paper we find that SO(10) with type II seesaw could be such a 
theory. 

When the fermion masses are given by the Eqs. (|2J) and (jl]), several possible outcomes are 
obtained by simple assumptions. To understand these possible outcomes from the Eqs.(j2]) 
and (jlj), let us first ignore h! . We then have the following possibilities: 

Assumption 1: 

Take h, f are general rank 3 matrices, and / is small. This is the case analyzed to fit 
observed experimental data and to obtain predictions from the minimality of the number of 
parameters in various papers 0]. Here, we list the properties resulting from the smallness 
of / without resorting to any numerical fit. 

• The CKM mixings are small, due to the fact that there is an approximate up-down 
symmetry and V C km = V u Vj 0. 

• Bottom-tau unification up to 0(f/h). 

• The 3 neutrino mixings are generically of 0(1) since h and / are unrelated matrices. 
The type II seesaw dominance of the neutrino mass is crucial for the generic largeness 
of the neutrino mixings. 

Thus it is interesting that without any special assumption, the gross features of fermion 
mixings can be reproduced. This does not, however, throw any light on the mass hierarchies 
among quarks and leptons, though one can fit the experimental results by the choices of 
parameters (even in the type I seesaw) 0, Q, [l5]. Since we use the experimental data as 
an input, these scenarios do not provide a fundamental understanding of either the mass 
hierarchy for quarks and charged leptons, or why the 13 neutrino mixing (U e z) is less than 
0(1). Similar situation holds for models where h! is added 
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Assumption 2: 

Let us next consider the specific case when h is a rank 1 matrix [7j, and / is a rank 3 
matrix with eigenvalues of / beinghierarchical (fi,f% <C i^) and small compared to the 



elects of H. As we noted in refft this cho.ce he.ps to supptess proton decay in SUSY 
SO (10) models without invoking huge cancellations among the colored Higgsino exchange 
amplitudes. In this case we will show in the next two subsections that the following results 
follow: 

1. CKM mixings are small. 

2. Approximate bottom-tau unification occurs. 
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3. — : — : — £ ~ r 2 : 1 : -3. 

m t rrib m T 

4. The quark mixing are related as V cb ~ m s /m b + e tcr m c /m t (where a is a phase) and 
V ub ~ Ki/j/Zs. 

5. Atmospheric and solar neutrino mixings are generically large, but 13 mixing is ~ fij f?>- 

All these predictions are in qualitative agreement with observations. The advantage of 
the rank one assumption is that it naturally explains the mass hierarchies among quarks and 
leptons in addition to large lepton mixings. We emphasize that these features are obtained 
from the rank 1 assumption above without using any numerical inputs, and our claim here 
is not based on the scenario of the numerical predictions from a fit in which the minimality 
of parameters plays a key role. 

Before we do a full demonstration of these results in the context of a three generation 
model, let us illustrate the first four points in the context of a two generation model. 



A. A two generation illustration 

In this subsection, we apply our rank one hypothesis to the second and the third gen- 
eration. We will confirm the results 1-4 mentioned above. The starting point is the mass 
relation from Eqs. ([2]) and (jl]) where we ignore the h' contribution . Using our assumption, 
we have h — ( sin 9 cos 9 )* ( sin 9 cos 9 ) h 3 and / = diag(/ 2 , / 3 ) (without loss of generality, 
we can parameterize / to be diagonal). The parameter 9 is of 0(1) in general. We now have 
ten parameters (9, h 3 and r± as real parameters, f 2 , f 3 , and r 2 as complex parameters, and 
vl for neutrino mass scale) describing ten observables of all lepton and quark mixings and 
masses. 

One can easily obtain r\m t ~ m b tdnf3 ~ m T tan/3 at the leading order neglecting 
0(f 3 /h 3 ) correction, where tan/3 is a ratio of vevs of H u and H d . Therefore, r x corresponds 
to the freedom of tan/3, t\ ~ tan/3/50. 

When / 2 <C /3 <C h 3 , we obtain 

m c h ■ 2 n m s h ■ 2 n m M oh ■ 2 a / K \ 
— ~ r 2 — sin 6, — ~ — sin 6, — - ~ — 3 — sin 9. (5) 

m t h 3 m b h 3 m T h 3 

Because m c /m t m s jm^ r 2 is small, i.e. r 2 ~ m c /m t / \m s /m b ). 

To proceed further, we first diagonalize the charged fermion mass matrices to zeroth order 
in f 2 , f 3 — > 0. The matrix diagonalizing this is given by 




tfo = . „ „ . ( 6 ) 
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since Uo(sm9 cos^) 4 = (0 1)*. Let us now see how the small quark mixings arise despite 
large mixings in Uq. Because UqYJJq has a small off-diagonal element, ri(f 2 — fa) sin 9 cos 9, 
the down-type quark mass matrix is diagonalized by V& = VdU , where Vd is close to unit 
matrix whose off-diagonal element is ~ f^/hs sin 9 cos 9. The up-type quark mass matrix is 
diagonalized by V u = V u Uo, where the off-diagonal element of V u is ~ r^fzjhz sin 9 cos 9. The 
quark mixing matrix is then given by Vckm = V u Vj = V u Vj , in this product U cancels out 
leaving small mixings between the two generations i.e. small V c b ~ {l — T^f^/h^ sin 9 cos 9 ~ 
(m s /mb + e lcr m c /m t ) cot 9, where a is a phase of r 2 . 

Coming now to lepton mixings, suppose that the charged lepton mass matrix is diag- 
onalized by Ve, then it can be written as Ve = V e Uo, where V e is close to a unit matrix 
similarly to the quark sector, and is roughly equal to V(t KM . Since the neutrino mass matrix 
is already diagonal as a parameterization, the PMNS matrix is given by the charged lepton 
mixings so that C/pmns = V e * — ^ckm^l- This leads to a large lepton mixing as desired. 
In the two generation case, 9 describes approximately (up to small corrections of order V c b) 
the "atmospheric mixing angle" . Since this was an input into our rank one ansatz, we can 
choose to be large to explain the observations. 

If /2; fs and r 2 are assumed to be real, there are six real parameters in this model. In this 
case, mb,m s , # atm and m 2 /m 3 can be written as a function of m c ,m t , m /J ,m r and V c b for 
example. Even if f 2 , fs and r 2 are all complex, we have the following approximate relation 
at the grand unified scale: 



Since the relations are satisfied under the assumption of approximate rank 1 property ir- 
respective of the counting of freedom, they are stable even in the case of three generation 
model. Indeed, the predictivity from the minimality of the parameter is related to the 
0(fs/h-s) corrections, and the minimality does not play a crucial role in the approximate 
relations from the rank 1 assumption. 

It is known that there is a solution that the large atmospheric mixing is obtained even if 
the smallness of / is not assumed a priori. In the scenario, the b-r mass convergence as well 
as the other experimental inputs predict the neutrino mixing as an output [4J. Our main 
goal in this section is not to give numerical predictions but rather to show how one can get 
qualitatively expected hierarchical pattern for masses and mixings. Later on we of course 
study the detailed numerical predictions. As it turns out there is a fine-tuned solution to 




(7) 
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fit the experimental data even if is comparable to /13. however such a fine-tuned solution 
is sensitive to the numerical inputs, and therefore the numerical predictions in this case 
may be unstable under a possible higher order correction. In our case where the quark and 
lepton mass hierarchy is predicted by the rank 1 assumption, they are stable under radiative 
corrections. 

While the qualitative predictions are in the expected range, we note that the approximate 
relation tan# atm ~ (fn s / m b)/Vcb is not very good agreement with the current observation, 
and small h! will be invoked to obtain the best fit of the experimental data. We emphasize 
that our final solutions do not use any fine tuned cancellations, and thus are stable even if 
we add small corrections to fit the numerical experimental data. 



B. Three generation case 



The fermion mass equations for this case are those in Eq. ([2]) with all coupling matrices 
being 3x3. The assumption that h has rank one means that we can write it as 



h 



b 

f = diag(/i,/ 2 ,/s 



c b a 



(/i,2 < /a)- 



(8) 



(9) 



Again, we can parameterize / to be diagonal and a, b, c to be real without loss of generality. 
At first, we ignore h! . In order to analyze the detailed consequences of this assumption, we 
go to the basis where h is diagonal. This is achieved by the matrix: 



( 



\ 



cos6 s sin9 s > 

cos 9 a sin 9 S cos 9 a cos 9 S — sin 9 a 
sin 9 a sin 9 S sin 9 a cos 9 S cos 9 a j 



(10) 



where tan 6^ = — c/b and tan# a = \/b 2 + c 2 /a with 

U hl4 = diag(0,0,/i 3 ), 



(11) 



where = a 2 + b 2 + c 2 . It is interesting to note that in the diagonalization matrix there is an 
ambiguity resulting from the residual SU(2) flavor symmetry in h (i.e. one of three mixing 
angles is not fixed at this stage). We choose the unitary matrix Uq to be an approximate 
leading order diagonalization matrix of Y u , Yd, and Y e as in the previous subsection (V u = 
V u Uq, Vd = VdUo, and Ve = V e Uo where V u ,Vd,V e are close to a unit matrix). Then, once 



S 



the / contribution is included, the afore mentioned SU(2) flavor symmetry is broken and 
the ambiguity in mixing angles alluded to above is removed. We wish to note that in our 
original parameterization of Uq, we chose the 13 element to be zero since even after including 
the /-contribution, the 13 element goes to zero in the limit — > (which is the limit 
where Y u ,d,e is rank 2). We have not used prejudices from neutrino experiment. 

By the same argument as in the case of two generations, Uq is cancelled out in the CKM 
mixing matrix and the quark mixings are small. The PMNS matrix is given by 

£/pmns = v;u , (12) 

and since the off-diagonal elements of Vi are small (being related to quark mixings), ne- 
glecting the 23 and 13 quark mixings, we get for the solar and atmospheric mixing angles 

E 

#atm ^ Oa, (13) 
#0 ~ 9 S ± #13 cot 9 a cos a, 

where a is defined as the diagonal phase matrix diag(l, e ia , e J/3 ) needed to diagonalize the 



charged lepton mass matrix 181 ] . 



We also get a formula for U e3 as follows: 

U e3 = (V e ) 12 sm9 a . (14) 

To proceed with the rest of the masses and mixings, let us define the matrices in the Uq 
rotation: Y a = UqYJJq, f = U^fU^ and so on. In this notation, V a is a diagonalization 
matrix of Y a . Because fa 3 = (fa — fa) sin 9 a cos 9 a and / 13 = (fa — f\) sin 9 a sin 9 S cos 9 S , one 
can obtain 

V ub ~ Vcb — • (15) 

J3 cos 9 a 

Neglecting 0(f 3 /h 3 ) and 0(fi t2 /h 3 ) corrections, V a can be approximately as: 

^ cos 9 — sin 9 
V = sin9 cos 9 , (16) 
{ l) 

where 

. x fa ~ fi cos 9 a sin 9 S cos 9 S 

sm0~ — -22 ' 17 

73 sm 9 a 



and thus 



fa ~ fi 

U e3 ~ — cot9 a sm9 s cos9 s . (18) 

73 
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Thus we have obtained all the features listed before. Due to the generic largeness of the 
relative mixing angles of the unrelated matrices, solar and atmospheric neutrino mixing 
angles are of 0(1) generically. On the other hand, 13 mixing is not in the category of the 
generic largeness since it is related to the ratio of eigenvalues of /. The eigenvalue ratio is 
also related to V u b/V c b implying that the 13 mixing angle has to be small in our approach. It 
is important to note that we do not assume a particular flavor texture such like hierarchical 
pattern in one matrix to obtain the feature. The key property to obtain the features for the 
neutrino mixings is that the correction to the rank one charged lepton mass matrix and the 
type II seesaw term are unified (or more roughly, simultaneously diagonalized) , as a result 
of SO(10) unification. 



C. Realistic model with h' 



The discussion above gives the qualitative consequences of the rank one property, and the 
experimental inputs are not used to obtain the features. The discussion below will address 
the issue of the experimental data for the first generation. Actually, we have not listed the 
first generation masses and V us before. In fact, if h! = 0, one obtains the following relation 
among the fermion masses: 

nhL : 1^ : UHH ~ r2 (l + r 2 X) : 1 + X : -3(1 - 3X), (19) 
m t nib m T 

where X = f\f2f^/{.o? fif'2 + b 2 f\f^ + c 2 f'2}^). When one fits down quark and electron masses 
(e.g., X = 0.35), the up quark mass is clearly too large since r 2 ~ m c /m t /(m s /mb) ~ 0.1. 
As a result, one of the first generation masses cannot be fitted. Besides, since V in Eq. ffT6]) 
is common for up- and down-type quarks, V us becomes too small compared to observations, 
since V us ~ V c bV u b. Therefore, one needs non-vanishing contribution from h! to obtain 
realistic masses for the first generation and V us under the rank 1 assumption. 



As is well-known, the empirical relation V us ~ \jTnd/m s is obtained when (Yrf)n — > and 
{Yd) 12 — {Yd) 21- Therefore we choose fn — > 0. When {Y d )n, {Y e )u — ^ is assumed, the 
choice of {Y e )i 2 {Y e )2i ~ {Yd) i2{Yd) 21 satisfies the Georgi-Jarskog (GJ) relation {m e m^m T ~ 
mdm s mb) for the down-type quarks and charged lepton masses. The up quark mass can be 
fit by using the freedom of r 3 . As a result, we have the following two solutions typically to 
fit the first generation masses and V us in a simple manner. 

Case A: f u ~ and |/ 12 + h' 12 \ — \ — 3/i2 + h' 12 \. The smallness of up quark mass is realized 
by a cancellation in {Y u )i 2 = r 2 f 12 + f^h' l2 . 

In this case, h! has to be symmetric which we can obtain by employing an extra 10 Higgs 
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field. For example, h' l2 — /12 is the simplest solution, giving 



U e 3 - -V us sm0 a , 

V us ~ 2sm0~2- — r^— -tanfc> s . 

/3 sin a 



(20) 
(21) 



where is given in Eq. (fT7|) and used a relation j\ ~ — / 2 tan 2 # s from / n ~ 0. Assuming that 
the corrections from the other elements of h! (e.g., h' 13 23 ) are small, we have an approximate 
relation 

(22) 



77^ ~ -V us tan 2 a , 

Vcb * 



which is in good agreement with the experiment. 

Case B: In this case we have fuf 22 — fu — 0. Then, 11 and 12 elements of {VfV t ) are zero, 
where V in Eq. (1161) is an approximate diagonalization matrix in the limit h! — > 0. The 12 
element of Vh'V 1 produces the Cabibbo angle. The GJ relation is manifest when \c e \ = 1 and 
(Vh'V^n — 0. The up quark mass is fitted by the smallness of r 3 , m u /m c ~ r 3 /r2md/m s . 

In this case, h! can be either symmetric or anti-symmetric. Since 11 element vanishes 
automatically, anti-symmetric coupling from 120 Higgs field is a better choice. As is noted, 
V contributes to U e3 , but it does not contribute to the Cabibbo angle. As a result, we obtain 
from Eq.flU]) 



\U e3 \ ~ sin0 a |sin0 + e n c e -K s | - 



f 2 1 

— tan S cot a + e n c e -K s sin a 

h 



(23) 



where 7 is a relative phase between / 12 and h' 12 roughly, and we have used a relation /1 ~ 
— ^tan 2 ^. Since V us is generated purely from h', it is not directly correlated to V uo /V c b 
contrary to the case A. 

As we have noted, to fit V c b, m s /mi, and atm very well, one needs a correction in h 23 . 
However, the correction does not affect the approximate expressions for U e3 very much. 

It is interesting that the U e3 is related to the mass ratio of neutrino in both cases. Since 
the GJ relation and the empirical relation of V us are not exact relations, there can be a shift 
from them in a numerical fit analysis. 

Here we assumed (V a )n — > to satisfy the GJ relation and the empirical relation of V us 
in a simple manner. When one introduces other parameters especially for symmetric h', 
there will be an accidental fine-tuned solution for the relations in a general fit for (y^ii 7^ 0. 
Actually, when h! is symmetric (c e = 1) and r 3 = 1, it results in a minimal model in which 



only one 10 and 126 Higgs fields couple to fermions with h being rank 3. In the minimal 
model for the fermion sector, it is known that there is a fine-tuned solution to fit fermion 



11 



masses and mixings 15[ unless the minimality of the Higgs potential is taken into account. 
In this case, when the first generation masses and V us are tuned, there is no freedom to 
adjust [7e3 and thus the approximate relation in the previous subsection holds, U e % ~ f2jfz- 
When r-i 7^ 1 (but r% ~ 1), U e s can be tuned to be any value (including zero) since first 
generation masses and V us can be fitted even if {Y e )i2 = 0. (When h is rank 1 and h! is 
anti-symmetric, there is no such fine-tuned solution. When h is rank 3, the fine-tuning fit 
for [7e3 = is allowed la].) Therefore, the assumption (Y a )n — > to satisfy GJ relation in a 



simple manner is crucial to keep the U e ^ prediction. Actually, when (Y a )n — ^ is assumed, 
the fine-tune solutions are removed, and the U e s is predicted as we have noted, irrespective 
of the number of parameters. We also note that the assumption (Y a )u — > is preferable to 
suppress nucleon decay amplitudes naturally. 

The case with the assumption that h (rank 1) and / are real and anti-symmetric h! is 
pure imaginary (in which case, the charged fermion mass matrices are hermitian) is in fact 
the model discussed in 0] . For this case, cancellation cannot happen between / and h' and 
thus the numerical fit does not shift very much from the above expression. In the numerical 
fit, it predicts \U e ^\ = 0.08 — 0.12 in the case where |c e | = 1 (and the GJ relation is manifest). 
Under the hermiticity assumption, one obtains e n = ±z and it is consistent with the above 
expression. In this case, since (Y e ) 2 2 is real, the PMNS phase is roughly same as the phase 
in the expression in Eq. (1231) and thus 

tantfpMNS - o— -r, ( 24 ) 
6 smu 

and we obtain <5pmns — ±30° or 180 ± 30° using the experimental inputs. 

For the case of |c e | ^ 1, one can also fit the experimental data and the prediction is |t/ e3 | = 
0.05 — 0.14. Since the cancellation is not allowed between / 13 and h ; 13 , the experimental data 
of V u b cuts the upper region of experimentally allowed mass squared ratio difference, and 
then gives an upper bound of U e3 . 



IV. TRI-BIMAXIMAL ANSATZ 



In the previous section, we incorporated large lepton mixings but their values were inputs 
into the theory. In this section, we consider special cases where the dominant part of the 
lepton mixing is in the tri-bimaximal form 19|. This would require special form for the rank 
one matrix h. We envisage that the rank one form for h as well as the matrix forms for / 
and h' come from some vacuum alignment of flavon fields, e.g., 20l. |21|. 

In the triplet flavon models, the 3x3 matrix can be expanded by the tensor products 
of the flavon fields when there are three independent flavons. The three flavon fields can be 
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expressed (without loss of generality, by making unitary transformations) as 

1 = (0,0,1), 2 = (O,a,6), fo = (c,d,e). (25) 
In general, there is no reason for the flavon vevs to be hierarchical, and the large neutrino 



mixings can originate from a ~ 6, c~<i~e 22J . The experimental result from the neutrino 



oscillation seems to imply a special alignment of flavon vevs rather than the generic largeness, 
namely 20], 

1 = (0,0,1), 2 = (0, -1, 1)/V2, 03 = (l,l,l)/v / 3. (26) 



The vacuum alignment can be obtained by imposing a discrete flavor symmetry 2lj, leading 
to tri-bimaximal neutrino mixings. 

It is worth pointing out at the beginning that the aligned flavon fields can be written in 
several ways by choice of coordinates (or by making unitary transformations) and the final 
results are independent of the coordinate choice. 

It is interesting to note that the aligned flavon vevs correspond to a link of a hexahedron, 
a diagonal line of a lateral surface, a diagonal line of a regular hexahedron, respectively. The 
interpretation becomes clear when the flavon fields are expanded in terms of the following 
orthogonal axes of coordinates (called hexahedral coordinate) 



xi = (1,0,0), x 2 = (0,1,0), x 3 = (0,0,1), 



(27) 



which correspond to the three lateral links of the regular hexahedron. The hexahedral 
coordinate is convenient to describe the Z 4 rotation around the surface-diagonal axes of 
the hexahedron. In fact, the regular hexahedron has different coordinates to describe the 
symmetry of the shape. One can consider a coordinate system which proves convenient to 
describe the Z 3 rotation around the vertex-diagonal axes of the hexahedron, 



(2,-l,-l)/V6, 



;i,i,i)/V3, 



(0,-1,1)7^2 



(28) 



The axes of the coordinates x\ and x' 3 are on the regular triangle which is formed by three of 
the hexahedron's vertices, and x' 2 is perpendicular to the triangle. We call this tetrahedral 
coordinate. 

The unitary matrix for the coordinate transformation from hexahedral (unprimed) to 
tetrahedral (primed) coordinate is the tri-bimaximal (TB) matrix i.e. x' 

( A fl ./I n \ 



xU^ B , where 



XJ- 



TB 



V 



(29) 



/ 



Therefore, in general (irrespective of the rank one assumption), if the charged- lepton mass 
matrix is (nearly) diagonal in the hexahedral coordinate and the neutrino mass matrix is 
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(nearly) diagonal in the tetrahedral coordinate, then the neutrino mixing matrix is (nearly) 
tri-bimaximal and given by: 

£4insp = VMtbV}, (30) 

where V e is a diagonalizing matrix of Y e in the hexahedral coordinate, and V u is a diagonal- 
izing matrix of Ai u in the tetrahedral coordinate. 

Suppose that the vacuum alignment in Eq. (l26|) is given in the hexahedral coordinate as 
we have noted. Then, those flavons in the tetrahedral coordinate are given as (0^ = 0i[/tb) 



[-1,V2, V3)/V6, 



(0,0,i; 



(0,1,0). 



(31) 



Therefore, from the discussion in the previous section, one can easily check that the nearly 
tri-bimaximal neutrino mixings are obtained when h is rank one formed by 0i (irrespective 
of the choice of the coordinate), and / is formed by 02 an d 03- We define 04 which is 
obtained an outer product of 2 and 3 , i.e. 4 = 3 x 2 . In the tetrahedral coordinate, 
4 = (1,0,0). 

As we have mentioned, the Yukawa matrices can be expressed in terms of the tensor 
products of the flavon fields. The symmetric matrices can be formed by six bases. Since we 
set 11 element to be zero in the hexahedral coordinate, we define the following matrices as 
the bases to form the linear space of symmetric matrices: 



Y 



>\9i 



^2 



/o o^ 
000 

v° l J 
/ 
1 

v° - 1 



\ 

-1 
1 

/ 



1 1 

1 i i 

22 

\ x 2 2 



\/6(0 



2^3 



+ 



^2) 



Y 5 = ^(0204 + 



U2) 



' 

-1 
1 

/ 
-1 

1 



/ 



-1 1\ 

-2 

2 / 

-1 1 \ 

1 
-I j 



(32) 



(33) 



(34) 



(35) 



(36) 
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where the elements of the matrices are presented in the hexahedral coordinate. 

In the following, we will consider two models: (I) V u is a unit matrix (/ is diagonal in 
the tetrahedral coordinate), (II) V v is close to a unit matrix. 



Model I : V„ 



The Model I can have both case A and case B solutions as described in the previous 
section. The case A solutions are, however, more natural in this model. In order to obtain 
such a solution, we employ additional 10 Higgs to obtain a correction matrix h', and hi is a 
symmetric matrix. 

We arrange the h, f, h' couplings as follows: 



h = h 3 Y x , 

f = h 3 e(Y 2 + XY 3 ), 

b! = h 3 e X p Y4 (or h 3 e X p Y 5 ) . 



(37) 
(38) 
(39) 



Then, since the ratio of eigenvalues of / is 1 : A : —A/2, we obtain Am^/Am^ = |A 2 . 
In the parameterization in the previous section, Uq in Eq. fllOl) is the tri-bimaximal matrix 
because a : b : c = : y/2 : — 1. 

The fermion Yukawa matrices are 



Y u = h + r 2 f + r 3 ti 
( 



e\(r 2 -r 3 p) eA(r 2 + r 3 p) 

eA(r 2 - r 3 p) r 2 e(l + |) + r 3 xe\p — er 2 (l — |) 

\eX(r 2 + r 3 p) -er 2 (l - f) 1 + r 2 e(l + |) - r 3 xeXp J 

Y d = nih + f + ti) 

( eA(l-p) eA(l + p) N 

eA(l - p) e(l + |) + xeXp -e(l - £ 



\ v eA(l+p) -e(l-f) 1 + e(l + f ) - xeAp y 
F e = n(/i-3/ + /i') 



/ 



rih 3 







eA(-3 - p) 



eA(-3+p) \ 



eA(-3-p) -3e(l + f) + xeAp 
v eA(-3 + p) 



3e(l 



36(1 - f ) 



(40) 
(41) 
(42) 
(43) 
(44) 
(45) 



1 - 3e(l + |) -xeXp ) 



where x = —2 when h! oc Y4, and x = 1 when /i' cx I5. 

For the numerical fits, the parameter e is given by e ~ m s /rrib ~ Vcfe, and the parameter 
A is given by A(l - p) ~ V us . 
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FIG. 1: f/ e 3 is shown as a function of Ara^/ Am^ trn (left) and # s (right) for Model I. 

When p ~ — 1, then Goergi-Jarskog relation is satisfied naturally. At that time, (Yd) 13 — 
0. This is interesting since the empirical relation \Vu\ — V us V c b is satisfied simultaneously 
when r 2 ,r 3 are small (up-type quark masses are more hierarchical rather than down- type 
quark masses). 

The parameter r 2 is fixed as |t~2 | ~ Tn c /mt/{m s /m^). The up quark mass can be made 
small by a choice r 3 ~ r 2 / p. 

Cabibbo angle, U e ^ and the ratio of mass squared differences Am^ ol / Am^ tm are all cor- 
related by the parameter A. The naive approximate relation is U e s ~ V us / (3\/2) as we have 
derived in the previous section. 

In Fig. 1, we plot U e % as a function of the mass squared difference ratio. In the plot, we 
fit m e /m T and m^/m T using p and e (which are assumed to be real in the plot). Then, U e % 
is calculated as a function of A. (The mixing angles do not depend on /i 3 ). We note that 
a correction from h' is needed to fit V c b, e.g., Ah' oc <j)\4>2 + 2 0i- As we have mentioned, 
such correction do not modify the U e z very much. In Fig. 1, we also show the plot of U e s as 
a function of 9 so \. Using the experimental constraint on Am^ ol / Aml tm and the other input 
from quark masses and mixings, we find that U e s is predicted to be 0.07-0.08. We note that 
the smaller side of experimental range of Am^ ol / Am1 tm is preferred from the numerical fit, 
which obeys from the naive relation A(l — p) ~ V us . The solar mixing angle 9 so \ is found to 
be ~ 32° from the plot which is obeyed by the approximation Eq. ( fTBl) when the PMNS phase 
is (or 7r), which is resulting from the assumption where p and e are real. The atmospheric 
mixing angle 8 atm is 45° up to ±2-3° correction from V c b for Model I irrespective of case A 
or case B solutions. 



The current allowed range for the neutrino parameters at 2a level are as follows 23]: 
#atm = 37° - 51°, 9 sol = 31.8° - 36.4° and Am 2 s J Aml tm = 0.027 - 0.038. 



1(3 



B. Model II: V v / 1 



The Model II can have both case A and case B solutions as well. In this model, the case 
B solutions are more natural. It is possible that the / coupling is not completely diagonal 
in the tetrahedral coordinate. Using the available freedom, we choose the 12 element of / 
in the hexahedral coordinate to obtain the case B solution. Since the 12 elements of Y 3 + y 4 
and Y3 + Y5 are zero, one can consider the choice: 



h = h 3 Y u 

f = h 3 e (Y 2 + X(Y 3 

h! = h 3 e\pYi. 



n)) 



(46) 
(47) 
(48) 



The matrix F 4 can be replaced with Y 5 . One can also choose h' to be antisymmetric, e.g., 

ti OC 2 3 - 030 2 . 

In the tetrahedral coordinate, the / coupling is written when / oc Y2 + X(Y 3 + Y4) (case 
Bl) as 



/ 



tetra 



V 










A 







5* 



\ 



/ 



(49) 



and if the notation in the previous section is used, we obtain a : b : c = (v^3cos^ + 
\^2smtp) : (\/2 cos tp — \/3 sin ip) : —1 where tan 2ip = y/E\/(l — A). The mass squared ratio 
is Am 2 sol /Aml tm = 3/4A 2 (l - 4A + 0(A 2 )). When we use / oc Y 2 + X(Y 3 + Y 5 ) (case B2), 



/ 



tetra 



-|A fx\ 



A 



f X 





1 



(50) 



/ 



we obtain a : b : c = (v^cos^' — sin^') : \/ 7 L : (—cosip' — \Z3sinip') where tan 2^' = 
v / 3A/(l + A/2). The mass squared ratio is Am 2 ol /Am 2 tm = 3/4A 2 (l - A + 0(A 2 )). 

We note that if there is a 0(A) correction in the 12 element in the tetrahedral coordinate, 
it modifies the 9 S angle largely, and it separates from the nearly tri-bimaximal mixing, and 
thus we do not use the choice. 

As we have obtained in the previous section, U e3 prediction is 



\U. 



e3| 



1 Am s 2 ol 



4m 



+ e 



3v/2 Ks 



(51) 



\| 2 Ami 

Clearly, if the parameters p, A, e are all real, then, 7 = or it, and the maximal and minimal 
values of U e3 are obtained. At that time, there is no phase in the PMNS mixing matrix. 
(The Kobayashi-Maskawa phase can be obtained from a phase of r 2 and/or r 3 .) 
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0.00 0.01 0.02 0.03 0.04 0.05 0.00 0.01 0.02 0.03 0.04 0.05 



Am 2 sol / Am 2 atm Am 2 sol / Am 2 atm 

FIG. 2: sin #i3 is shown as a function of Am^/Am 2 ^ (left) and # a t m is shown as a function of 
^ m soi/^ m atm (right) for Model II-case Bl (described in the text). 




Am 2 sol / Am 2 atm e atm (degree) 

FIG. 3: #atm is shown as a function of Am^/Am^ (left) and sin 6*13 is shown as a function of 
#atm for Model II-case B2 (described in the text). 

In Fig. 2 (case Bl, Eq. (l4"9l) ) and Fig. 3 (case B2, Eq. (l50]) ). we plot U e ^ when p A, e are real 
to find the lower and upper limits. Due to the off-diagonal elements of / in the tetrahedral 
coordinate, the atmospheric angle shifts from 45°, and the shift is correlated to U e 3, unlike 
the case of Model I. 

In case Bl, using the experimental constraint on Am^/Am^, we find that U e ^ is 
predicted to be 0.05-0.08. This solution corresponds to the sign choice e n = —1 in Eq. (j5"TT) . 
We also plot atmospheric mixing angle in Fig.2 (right). It is interesting to note that for the 
case Bl, A should be negative for |A| ~ 0.1 — 0.3 to fit mass squared difference ratio since 
^ m soi/^ m atm = 3/4A 2 (l — 4A + 0(A 2 )). As a result, the direction of the shift is determined 
to fit experimental values, i.e. # atm > 45°. 
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35.0 




Am 2 sol / Am 2 atm 




0.05 



Am 2 S0 , / Am 2 



atm 



FIG. 4: 9 so \ is shown as a function of Am^|/Amj tm for Model II-case Bl (left), case B2 (right). 



In case B2, the mass squared ratio can be fitted for both signatures of A (we find two 
branches in both graphs in Fig.3). The plot in Fig. 3 is shown constraining tan 2 8 so \ > 0.35, 
and one can find that Am^/Am^ becomes too small for one of the branch to fit solar 
mixing angle, and # atm > 45° is favored in this CclSG clS well. Larger U e3 values (> 0.15) is 
preferred once we include the experimental limit on Am;? ol /Am 2 tm . The solution corresponds 
to the sign choice e* 7 = +1 in Eq. (151jl . 

As we discussed above that both case A and case B solutions can be obtained in Models 
I and II. The interesting question is how can we distinguish these two models. For example, 
the predictions of U e % (as shown in Figs 1,2 and 3) distinguishes between the cases A and B. 
Actually, if U e % is just below the current CHOOZ bound, the case B solution with e* 7 = 1 
is preferred. Since the case A is a more natural solution for Model I and case B is a more 
natural solution for Model II, one may weakly conclude that this prediction distinguishes 
between Models I and II. But a stronger way to distinguish these models would be to use 
the predictions of 9 atm . In Model I, 8 at m is fixed to be 45° up to ±2 — 3° corrections from V c b, 
where as Model II prefers 8 atm > 45° (Figs 2 and 3) in the experimentally allowed region. 
This difference is directly due to the rigorousness of fi — r symmetry in the / coupling from 
the tri-bimaximal ansatz (V„ = 1 (Model I)) where the / coupling has Z 2 x Z 2 symmetry. 
In Model I, the deviation from the maximal angle is related to V c b, while in Model II, the 
deviation is related to U e s. The current best fit value of the atmospheric mixing angle 
(6»^ tfit = 43° Q) is nearly the maximal mixing, and it implies the Model I. However, the 
error is still large. The accurate deviation from the maximal angle will be obtained in future 
three generation fit of the neutrino oscillations 24(, and it will give us an important test for 
the tri-bimaximal ansatz. 

The predictions for 6 so \ are similar with a small margin in these models, since the devia- 
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tions from the tri-bimaximal angle (8 S = 35.3°) are related to U e ^ in all cases. In Model I, 
# so i is predicted to be ~ 32°, where as, in Model II (Fig. 4), ao \ is predicted to be ~ 34° 
(case Bl), ~ 33° (case B2), once we include all the experimental bounds. The PMNS phase 
is assumed to be or ir in the plot, and the general phase fit will change the predictions of 
the angle, especially for Model II. 



V. DERIVATION OF RANK ONE ANSATZ 

The rank one Yukawa coupling with 10 Higgs field generates the features of flavor hi- 
erarchy, and rank 1 matrices can often appear in various ways (flavor symmetry, discrete 
symmetry, and string models). In this section, we give an SO(10) model, where the rank 
one ansatz used in our discussion of flavor emerges from a discrete symmetry. 

When the direct couplings of chiral fermions with a Higgs field are forbidden by a sym- 
metry, and the effective Yukawa couplings are generated by propagating vector-like matter 
fields, the rank of the effective Yukawa matrix depends on the number of the vector-like 
fields. Actually, when there are only one pair of vector-like matter fields as a flavor singlet, 
the effective Yukawa matrix is rank 1. 

The model we assume has one extra vector-like pair of matter fields with mass slightly 
above the GUT scale (denoted by ijjy = 16y © z/v = 16y) an d three gauge singlet fields 
Y a . We add a Z4 discrete symmetry to the model under which the fields if) a —>■ iip a , and 
Y a — > —iY a . The 10-Higgs field H is invariant under this symmetry. The gauge invariant 
Yukawa superpotential under this assumption is given by 

W = ^vHX^Jv +M v VvVv + Vv£Y a ?/v (52) 

a 

When we give vevs (Y a ) ^ 0, ipy and ip a are mixed. The heavy vector-like fields, ipy and a 
linear combination of Vv an d ip a (he. Myipy + J^a^a^a), and the effective operator below 
its scale and at the GUT scale is given by: 



r x 



My + E a Y a 2 



H 



(53) 



This gives rise to a rank one h coupling. We note that it does not contradict the 0(1) top 
Yukawa coupling, when M v ~ J2 a Y a ( or M v < T,a Y a)- 

If we let the 126 Higgs field transform like —1 under Z 4 , it can induce the / coupling 
with rank three. 

Another way to get mass matrix patterns of types in sec. IV is to assume that there 
are three component flavon fields (denoted by the dimensionless field 0j = J|) which are 
representations of some internal flavor group and constrain their couplings to fermions by 
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other symmetries. As an example, let us choose four flavon fields which transform as follows 
under a Z 4 x Z 2 x Z 2 group: (f>i(-i, +, +); <f> 2 (l, -, +); 3 (1, +, -); 04(1, -, -); Hh +, +); 
A(— 1,+, +); H(l,+,+); H'{— 1, — , — ). The invariant Yukawa coupling under these sym- 
metries is: 



+ (fofaWA + 3 3 ^A + (fcfaWA + (fefaipipH' + h.c. (54) 



In general the vevs of the 0; fields align as in Eq. (25) which can then lead to our type of rank 



one models. With suitable discrete symmetries, e.g., A (27) j2lj, the flavon vevs can align 
as in Eq. (26), leading to models of the type considered here (Model I). When Z 2 group is 
chosen instead of Z 2 x Z 2 group such that 0i(+); (f>2(—)', 03(— ); 04(+); ^(+)> A(+); H(+); 
H'(+), 0203 term is allowed in the ipipA coupling and Model II (Bl) can be considered. 
More details on the flavon vev alignment with discrete symmetries and implications for rank 
one models is currently under investigation. 

VI. CONCLUSION 

In conclusion, we have shown how a simple ansatz for the dominant Yukawa coupling 
matrix in renormalizable SO(10) models can lead to a unified understanding of the diverse 
quark and lepton flavor hierarchies. We suggest this as a possible way to address the 
challenge of a unified description of quark-lepton flavor. We have not attempted in this note 
to derive our ansatz from any specific discrete symmetries, although we show a guideline to 
obtain a rank one form from a higher than GUT scale theory. This may be the next step 
towards a complete theory of flavor. 

Within our rank one hypothesis, we have considered two classes of models which are nearly 
tri-bimaximal and point out that both a measurement of # 13 and the atmospheric mixing 
angle 8 atm can distinguish between these models. In these models, the natural minimal value 
of #i3 is around 0.05, whereas the maximal value can be larger than 0.15. This range of #13 
can be probed in the upcoming experiments. The current estimate of #13 using the 1.5 a 



excess of events in MINOS — v e appearance channel 25j and all other experimental data 
is sin 2 13 ~ 0.02 ± 0.01 (1 a) [26]. 
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